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ABSTRACT 
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Aims. We report the discovery of CXO J1415. 2+3610, a distant (z ~1.5) galaxy cluster serendipitously detected as an extended source 
with a very high significance level (S/N ~ 1 1) in a deep, high-resolution Chandra observation targeted to study the cluster WARP 
J1415. 1+3612 at z = 1.03. This is the highest-z cluster discovered with Chandra so far. Moreover, the total exposure time of 280 ks 
with ACIS-S provides the deepest X-ray observation currently achieved on a cluster at z > 1.5. 

Methods. We perform an X-ray spectral fit of the extended emission of the intracluster medium (ICM) with XSPEC assuming a 
■ single-temperature thermal mekal model. We use optical and infrared (IR) observations from Subaru-5;(/>n'w?e (BVRiz), MOIRCS 

(JKs), and Spitzer/IRAC (3.6 /jm) to confirm the presence of an overdensity of red galaxies matching the X-ray extended emission. 
i-C ' We use optical and IR data to investigate the color-magnitude relation of the candidate member galaxies. 

Oh. Results. From a preliminary X-ray spectral analysis, we detect at a 99.5% confidence level the rest frame 6.7-6.9 keV Iron K a line 

complex, from which we obtain zx = 1-46 ± 0.025. Our X-ray redshift measurement is supported by the optical and IR data. The 
analysis of the z - 3.6//m color-magnitude diagram shows a well-defined sequence of red galaxies within 1' from the cluster X-ray 
emission peak with a color range [5 < z-3.6/im < 6]. The photometric redshift obtained by spectral energy distribution (SED) fitting is 
Zphot = 1.52 + 0.06. After fixing the redshift to z = 1.46, we perform the final spectral analysis and measure the average gas temperature 
with a 20% error, kT = 5.8*}'? keV, and the Fe abundance Z Fe = 1.3*n!Z . We fit the background-subtracted surface brightness with 
a single beta-model out to 35 arcsec (the maximum radius where the X-ray emission is detected), and derive the deprojected electron 
density profile. The ICM mass is 1.09*2| x 10 13 M G within 300 kpc. Under the assumption of hydrostatic equilibrium, the total mass 
is M 25 oo = 8.6+ 2 ^ x 10 i3 M Q for R 250 o = (220 ± 55) kpc. Extrapolating the profile at larger radii, we find M S0Q = 2.1*ol x 10 14 M Q 
for R S00 = 510!;^ kpc. This analysis establishes CXOJ1415.2+3610 as one of the best characterized distant galaxy clusters based on 
X-ray data alone. 

in : 

^s^j i Key words. Galaxies: clusters: intracluster medium; individual : CXO J1415. 2+3610 - X-ray: galaxies: clusters 
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1 . Introduction of the baryonic fraction as a function of the cosmic epoch and 

therefore does not depend on the density evolution of galaxy 

Clusters of galaxies are both important probes to measure the clusters. 

growth of cosmic structures of the Universe and excellent as- There is now a general convergence towards the standard 
trophy sical laboratories (see Kravtsov & Borgani 2012, and ref- A cold dark matter model (ACDM), which is based on ob- 
H 1 erences therein). Their measured number density and evolution servables related to the cosmic microwave background (CMB, 
. . .' strongly depend on the cosmological parameters, particularly at Komatsu et al. 201 1), the correlation properties and the growth 
the high mass end (see Allen et al. 2011). Their number den- rate of the large-scale structure of the Universe (Reid et al. 
sity evolution on a time scale of ~ 9 Gyr, corresponding to 2012), and the supernova (SNae) surveys (for recent develop- 
the highest redshift where virialized clusters are currently found ments see Suzuki et al. 2012). Although CMB and SNae ex- 
(z > 1.3), ranges between one and more than three orders of periments that directly probe the geometry of the Universe 
magnitudes for masses M > 5 x 1O 14 /j _1 M , depending on the have gained large visibility, significant efforts are still being de- 
cosmological parameters (see Rosati et al. 2002). Therefore, the voted to cosmological tests with galaxy clusters thanks to the 
density of massive clusters at high redshift allows one to trace high complementarity with the other tests (Mantz et al. 2010). 
the growth rate of perturbations, thus providing a robust con- Recently, the emphasis has shifted towards the search of devia- 
straint for cosmological models. Other cosmological tests with tions from the reference ACDM model, like tests of general rela- 
galaxy clusters are provided by the power spectrum of galaxy tivity (Schmidt et al. 2009; Rapetti et al. 2012), coupled dark en- 
clusters (Balaguera-Antolmez et al. 201 1) and the baryonic frac- ergy models (Amendola et al. 2012), or non-Gaussianity in the 
tion (Allen et al. 2002; Ettori et al. 2003, 2009). The latter is a primordial density perturbation field (e.g., Sartoris et al. 2010; 
geometrical test that relies on the assumption of the constancy Sawangwit et al. 2012). Potentially, the mere presence of few 
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massive, high-z galaxy clusters such as those recently discovered 
(e.g., Jee et al. 2009; Foley et al. 201 1; Santos et al. 201 1) could 
create some tension with the standard ACDM and quintessence 
models (e.g., Jimenez & Verde 2009; Baldi & Pettorino 2011; 
Chongchitnan & Silk 2012; Mortonson et al. 201 1; Hoyle et al. 
201 1; Harrison & Coles 2012; Waizmann et al. 2012). So far, no 
evidence for such a discrepancy has been found, except when 
considering the combined probability of the most massive, high- 
z clusters (Jee et al. 201 1). 

X-ray selected clusters play a key role in this context. Thanks 
to the X-ray thermal emission from the intracluster medium 
(ICM), which is the largest baryonic component of galaxy clus- 
ters, it is possible to measure the cluster mass with good accu- 
racy up to high-z. Unfortunately, X-ray studies of distant galaxy 
clusters are limited not only by relatively small samples, but also 
by low fluxes and surface brightness, which hamper a meaning- 
ful X-ray spectral analysis in most cases. This can be achieved 
only with very deep exposures, possibly with the high imaging 
quality that only Chandra can reach. An independent strategy is 
to find clusters and estimate their masses through the Sunyaev- 
Zeldovich (SZ) effect. Surveys carried out with the Atacama 
Cosmology Telescope (Sehgal et al. 2011) and the South Pole 
Telescope (Reichardt et al. 2012) are delivering the first remark- 
able results. However, the calibration of mass proxies for SZ- 
selected clusters still relies on a cross-calibration with X-ray 
data (Andersson et al. 2011; Planck Collaboration et al. 2012). 
Before SZ mass measurements are robustly calibrated, the two 
major existing X-ray facilities, XMM-Newton and Chandra, still 
offer the best way to find and characterize distant massive clus- 
ters. 

Significant progress in this field can be pursued with a two- 
step strategy. The first step consists of exploiting the Chandra 
and XMM-Newton archives to identify high-z cluster candi- 
dates. The total solid angle covered by both telescopes amounts 
to about 100-200 deg 2 (excluding the Galactic plane), and only 
a minor fraction of these fields reaches the required sensitivity 
to detect clusters at z > 1, whose typical flux is about 10~ 14 erg 
s cirT 2 . This has been best achieved by a number of ongo- 
ing surveys, namely, the XMM Distant Clusters Project (XDCP, 
Fassbenderet al. 2011a), by far the most sucessful to date, the 
XMM-Cluster Survey (XCS, Lloyd-Davies et al. 201 1), and the 
XMM Large Scale Structure Survey (XMM-LSS, Adami et al. 

201 1) 1 . As a second step, it is necessary to perform in-depth 
studies of the most massive, distant cluster candidates. This strat- 
egy requires time-expensive observations with Chandra in order 
to obtain good signal-to-noise ratio (S/N) spectra (typically col- 
lecting a total of 1000 net photons or more in the 0.5-7 keV band) 
and high-resolution images to remove the effect of contaminat- 
ing active galactic nuclei (AGN) emission and of central cool 
cores. Clearly, this goal can also be pursued with deep X-ray 
follow-up of SZ-selected targets. 

If we focus on the most distant cluster candidates, we find 
that only nine clusters have been spectroscopically confirmed at 
z > 1 .5 to date and only some of them have estimated masses in 
excess of 10 14 Mq (Brodwin et al. 201 1, 2012; Fassbender et al. 
2011b; Nastasietal. 2011; Santos et al. 2011; Papovich et al. 
2010; Gobatetal. 2011; Zeimannetal. 2012; Stanford et al. 

2012) . We note that most of them have been identified in in- 
frared (IR) surveys or by combining near infrared (NIR) and X- 



1 For a quick overview of current X-ray surveys, see Tundo et al. 
(2012) 



ray data 2 . Moreover, the apparent lack of extended X-ray emis- 
sion in some of these systems suggests low X-ray luminosities 
and points towards a re-classification of these objects as proto- 
clusters. Typically, data used for detection do not allow one to 
determine the virialization status of a structure, or to measure 
the average temperature of the ICM. 

In this paper, we present the identification and characteriza- 
tion of the most distant X-ray cluster discovered with Chandra 
at z ~ 1.5 in the same field of the cluster WARPJ1415. 1+3612 
(hereafter WARPJ 14 15, at z spec = 1-03, see Santos et al. 2012). 
Since these X-ray data are the deepest observation ever car- 
ried out with Chandra on clusters at z > 1 (370 ksec), they 
allowed us not only to firmly detect the extended ICM emis- 
sion, but also to measure its properties from the X-ray spec- 
trum. In addition to the X-ray data, we have high-quality multi- 
band optical/NIR data from Subaru and infrared Spitzer-IRAC 
(3.6fim) imaging. Currently, the only comparable systems in 
terms of redshift and estimated mass are XMMUJ2235. 3-2557 
at z = 1.393 (Rosatietal. 2009), XMMUJ0044 at z = 1.578 
(Santos etal. 2011), and, recently, IDCS J1426+3508 at z = 1.75 
(Stanford et al. 2012). Only the first one, however, was observed 
with sufficiently deep Chandra and XMM-Newton data to derive 
a robust total mass estimate. 

The paper is organized as follows. In Section 2 we describe 
our X-ray spectral analysis, which provides us with X-ray red- 
shift, temperature, and metallicity. We also measure the surface 
brightness profile and concentration and derive the ICM mass 
and the total cluster mass. In Section 3 we describe the optical 
and IR data, discuss the red sequence, and compute galaxy pho- 
tometric redshifts. In Section 4 we discuss the properties of CXO 
J1415.2+3610in a cosmological context and in comparison with 
other distant X-ray clusters. Finally, our conclusions are summa- 
rized in Section 5. Throughout the paper, we adopt the seven- 
year WMAP cosmology (Q A = 0.73, Q m = 0.27, and H = 70.4 
km s _1 Mpc"" 1 (Komatsuet al. 2011). In this cosmology, 1" on 
the sky corresponds to about 8.56 kpc at z = 1.5. Quoted errors 
always correspond to a 1 cr confidence level. Filter magnitudes 
are presented in the AB system, unless stated otherwise. 



2. Cluster detection and X-ray data analysis 

2.1. X-ray data 

CXO J1415.2+3610 (hereafter CX01415, RA: 14:15:20.143, 
DEC:+36: 10:27.56) was serendipitously detected in a deep 
Chandra GO observation of 280 ks with ACIS-S granted in 
Cycle 12 (PI J Santos) to study the ICM properties of the cool- 
core cluster WARPJ 14 15 at z — 1.03. The new system, lo- 
cated at a distance of ~ 2 arcmin in the south-west direction 
of WARPJ1415, was immediately detected by visual inspection 
as an extended X-ray source in the ACIS-S image (see Figure 1). 
The cluster was also visible in a shallower ACIS-I image of 90 
ks available in the Chandra archive, which we also use in this 
study. We tried to recover the XMM data for WARPJ 14 15 to ex- 
tract more information on CX01415. Unfortunately, CXO 14 15 
is barely visible in the total XMM image (PN plus 2 MOS) and is 
partially affected by a gap. The XMM exposure time after clean- 
ing from high-background interval is about 16 ks. Therefore, 
XMM data do not provide any useful signal, and our analysis 
is entirely based on the Chandra data. 



2 No SZ-selected clusters are currently found above z = 1.5, the 
most distant confirmed SZ-selected cluster is reported at z - 1.32 
(Stalder et al. 2012). 
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Fig. 1: ACIS-S+I image of the clusters field, smoothed with a 
3" Gaussian kernel. The bright cluster WARPJ1415 at z = 1 is 
visible in the upper right, whereas CX01415 lies at ~ 2' in the 
south-west direction. The solid circles have a radius of 35" (cor- 
responding to 300 kpc at z = 1-46), whereas the dashed circle 
with r = 24" refers to the aperture of maximum S/N used for the 
spectral analysis. The image has a size of 4' x 4'. 
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Fig. 2: Cstat as a function of redshift parameter for the combined 
(ACIS-S+ACIS-I) spectral analysis of CX01415. The minimum 
is obtained for z = 1 -46, when AC stat ~ 8. The cyan stripe shows 
the 1 cr confidence interval obtained from the photo-z (z p hot = 
1.52 + 0.06, see Section 3.2) 



The total exposure time of the analyzed data is 370 ks 
(ACSI-S+ACIS-I). The data reduction and source extraction 
were performed with CIAO 4.3 and CALDB v4.4.5. A detailed 
description of the data reduction procedures can be found in 
Santos et al. (2012). The circular region of maximum S/N is ob- 
tained for a radius of 24 arcsec centered on RA=14: 15:20. 143, 
DEC=+36: 10:27.56. The source is detected with a S/N of 11.6 
and ~ 3, corresponding to about 540 and 115 net counts (0.5-7 
keV band) in ACIS-S and ACIS-I, respectively. 

2.2. Spectral analysis: temperature, metallicity, and redshift 

We perform a combined spectral analysis of the ACIS-S and 
ACIS-I data with Xspec vl2.6 Arnaud (1996). The adopted 
spectral model is a single-temperature mekal model (Kaastra 
1992; Liedahl et al. 1995), using as a reference the solar abun- 
dance of Asplund et al. (2005). The local absorption is fixed to 
the Galactic neutral hydrogen column density measured at the 
cluster position and equal to Nh = 1.05 x 10 20 cirT 2 taken from 
the LAB Survey of Galactic HI (Kalberla et al. 2005). The fits 
are performed over the energy range 0.5 - 8.0 keV. We used 
Cash statistics applied to the source plus background, which is 
preferable for low S/N spectra (Nousek & Shue 1989). 

First, we search for the K B Fe line complex in order to 
measure the redshift zx directly from the X-ray analysis. For 
Chandra observations, about 1000 net counts are needed to mea- 
sure zx at a 3<r confidence level, and even more are needed in the 
case of hot clusters, as shown in Yu et al. (201 1). Nevertheless, 
we run a combined fit on the ACIS-S and ACIS-I spectra, leaving 
free redshift, temperature, metal abundance, and normalization. 
We find that the best fit gives zx — 1 -46 + 0.025. To better quan- 
tify the robustness of the detection of the Fe K a line complex, 



we plot in Figure 2 the quantity C stat as a function of the redshift 
parameter, while all the other spectral parameters are set to the 
best-fit values. The best-fit redshift is obtained with a AC stat ~ 8. 
Formally, a 3 cr confidence level is obtained when AC sla , = 9 
(see Yu et al. 201 1). As a result, the redshift measured from the 
X-ray analysis corresponds to a ~ 2.8<x (99.5% confidence level) 
detection of the Fe line. However, as we show later, another rel- 
evant piece of information comes from the photometric redshift, 
which allows us to exclude secondary minima at z ~ 0.8 and 
z ~ 2.5 and increases the robustness of the redshift value ob- 
tained from the X-ray analysis. Therefore, in the following anal- 
ysis we assume z - Zx - 1-46. This result shows that the mea- 
sure of the redshift from X-ray spectral analysis is possible up to 
the highest redshift where X-ray clusters are currently detected. 
This aspect is extremely important when designing a wide and 
deep survey mission for the future, such as the Wide Field X- 
ray Telescope (Murray & WFXT Team 2010; Rosati et al. 2010; 
Tozzi et al. 2010). 

The cluster unabsorbed soft-band flux within a circular re- 
gion of 24" radius is 5o.5_2.oteV = (6.4 + 0.8) x 10~ 15 erg s -1 
crrT 2 , and the hard flux SViofev = (5.8 + 0.7) x 10~ 15 . The 
cluster emission-weighted temperature is kT = 5.8^[ q keV. The 
measured Fe abundance, in units of Asplund et al. (2005), is 
Zf e = 1.3^QgZ/r e0 . We warn that Fe abundance values measured 
freezing the redshift parameter to the zx value are slightly biased 
high (see Yu et al. 201 1). Finally, we note that this analysis is in 
very good agreement with that based on the ACIS-S data only, 
which indeed carries most of the signal. The inclusion of the 
ACIS-I data helps to slightly reduce the statistical error bars. For 
sake of clarity, we show in Figure 3 the spectrum of CX01415 
obtained only with ACIS-S, along with the best-fit model. 
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CX01415, z_X=1.46 




Energy (keV) 

Fig. 3: Folded Chandra ACIS-S spectrum of CX01415 with its 
best-fit mekal model (solid line). The iron K a line complex is 
detected at the 99.5% confidence level at 2.8 keV in the observed 
reference frame. 



2.3. Surface brightness profile and concentration 

To measure the surface brightness of CX01415 we use only 
the ACIS-S image. The azymuthally averaged surface brightness 
profile, which was computed in the exposure-corrected ACIS-S 
image in the soft band, is traced out to about 35" and is well de- 
scribed by a single B model, S(r) = S (l + (r/r c ) 2 ) _3 ^ +0 - 5 + bkg 
(Cavaliere & Fusco-Femiano 1976). The best-fit model, shown 
in Figure 4, left-hand panel, has a slope B - 1.2 and a core ra- 
dius r c = 240 kpc for a reduced x 2 — 0.7. The uncertainties 
on the profile parameters are quite large, given the degeneracy 
between B and r c (see Figure 4, right-hand panel). This degen- 
eracy does not significantly affect the overall profile within the 
extraction radius, which is used to derive both total and baryonic 
mass within 300 kpc in Section 2.4. However, the extrapolation 
of the surface brightness profile at larger radii is uncertain, and 
further assumptions are needed to estimate the total mass up to 
the virial radius. The total luminosity, being dominated by the 
inner regions, will be mildly affected. 

The luminosity emitted within 24 arcsec in the soft band is 
L x [0.5 - 2.0] = (0.61 + 0.07) x 10 44 erg s" 1 and the bolometric 
luminosity Lj,„i = (2.1 ± 0.2) x 10 44 erg s _1 . To compute the 
total cluster rest-frame luminosity, we apply a correction factor 
obtained by extrapolating the fit of the surface brightness profile 
up to the virial radius (see Section 2.4). We obtain a total soft- 
band luminosity of L x [0.5 - 2.0] = (0.93 + 0.10) x 10 44 erg 
s _1 and a bolometric Lt, i = (3.2 + 0.3) x 10 44 erg s _1 . This is 
probably an upper limit, given that surface brightness profiles are 
generally observed to steepen in the outer regions with respect to 
the extrapolation of the beta-fit model (Ettori & Balestra 2009). 

The rather flat surface brightness profile within 100 kpc does 
not support the presence of a cool core. To quantify the core 
strength, we compute the value of csb, a concentration parameter 
defined in Santos et al. (2008) as the ratio of the surface bright- 
ness within 40 kpc and 400 kpc: c SB - SB (< 40kpc) / SB(< 
4-OOkpc). A cluster is expected to host a cool core if c$b > 0.075. 
This simple phenomenological parameter has proven to be a ro- 
bust cool core estimator, particularly in low S/N data typical of 
distant clusters. We measured c$b — 0.059+0.017, thus confirm- 
ing that CX01415 does not show evidence of a well developed 
cool core. 



2.4. Mass measurement 

The baryonic mass can be directly computed once the three- 
dimensional electron density profile is known. For a simple beta 
model, this is given by n e (r) = n e o[(l +(r/r e ) 2 ) -3 ^ 2 ]. To measure 
n e o, we use the relation between the normalization of the X-ray 
spectrum and the electron and proton density in the ICM for the 
mekal model: 



Norm 



10 



4ttD 2 (1 +z) 2 



nndV , 



(1) 



where D a is the angular diameter distance to the source (cm), 
and n e and (cm -3 ) are the electron and hydrogen densities, 
respectively, and the volume integral is performed over the pro- 
jected region used for the spectral fit (24 arcsec). For CX01415, 
this gives a central density of n e o = (0.99 + 0.13) x 10~ 2 crrT 3 . 

The total cluster mass can be estimated under the assump- 
tions of hydrostatic equilibrium and spherical symmetry, which 
leads to the simple equation (Sarazin 1988) 

M{r) = -4.0 x 10«M e T(keV)r(Mpc{ + ^] . 

\ dlogr dlogr ) 

(2) 

Since we are not able to measure the temperature profile, 
we assume isothermality inside the extraction radius and be- 
yond that adopt a mildly decreasing temperature profile kT oc 
r~°- 24 , as found in local clusters (see Leccardi & Molendi 2008). 
Therefore the term in parenthesis in Equation (2) is contributed 
by a term due to the slowly declining temperature profile and by 
the slope of the density profile, which is simply -3Bx 2 /(l +x 2 ), 
where x = r/r c . Since we can robustly measure the slope of the 
density profile only up to 300 kpc, we decide to keep the slope 
constant beyond 300 kpc (dlog(n e )/dlogr = 2.2). 

To measure the total mass at a given density contrast, we 
solve the equation M^(r^) = 4/37ir 3 Ap c (z c i) to compute the ra- 
dius where the average density level with respect to the critical 
density p c (z c i) is A. Typically, mass measurements are reported 
for A = 2500, 500, and 200. The lcr confidence intervals on the 
mass are computed by including the error on the temperature 
and on the gas density profile. The results are shown in Table 
1. Only the radius r25oo = (220 + 55) kpc is well within the 
detection region. We obtain M 25 oo = 0.86+ ,; 2 } x 1O 14 M . We 
also report the total mass within 300 kpc M(r < 300kpc) = 
1. 38*0 23 x 10 14 M o and the corresponding ICM mass M ICM {r < 
300kpc) = 1.09!°; 3 |J x !0 I3 M o . 

The total mass extrapolated to r$oo, with the additional as- 
sumption about the temperature and density profiles, are M500 = 
2.07!° ™ x 1O I4 M for R 500 = 510+f kpc. We also extrapolate 
the mass measurement up to the nominal virial radius, finding 
M200 = 2.8+$ x 10 14 M o for ^200 = 760^° kpc. The ICM mass 
fraction is almost constant with a value ficu — 0.10 + 0.02, 
which is in very good agreement with the empirical relation be- 
tween ficM and M500 shown in Vikhlinin et al. (2009). 

Finally, for completeness, we recompute the masses after re- 
laxing the assumption of isothermality within the extraction ra- 
dius and adopting the temperature profile kT oc r ~ 0M at any 
radius. In this case, the only requirement is that the average tem- 
perature within this radius is 5.8 keV as observed. With these 
assumptions, masses at M500 and M200 are ~ 18% lower than 
those obtained with the simplest choice T = const for r < r exl . 
This difference should be regarded as systematic, since the tem- 
perature profile within the extraction radius may well behave dif- 
ferently from a simple power law, with significantly different ef- 
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Fig. 4: Left: Background-subtracted surface brightness profile in the soft (0.5-2 keV) band (points) and best-fit /3 model (red dashed 
line) for CX01415. Error bars correspond to 1 <x uncertainty. Right: 1-2-3 <x confidence levels for two relevant parameters in the 
f3-r c parameter space. 



Table 1: Summary of the mass estimates for CX01415. The "*" 
indicates that the mass values are extrapolated beyond 300 kpc, 
the maximum radius where the ICM emission is actually de- 
tected. 



Table 2: Total mass at R^oo and R200 obtained under dif- 
ferent assumptions for CX01415. (2) Vikhlinin et al. (2009); 
(3a) Vikhlinin et al. (2009); (3b) Fabjan et al. (2011); (4) 
Navarro et al. (1996). 



Radius M lol M IC m Method R soo M soo M 2( 

kpc M Q M Q kpc M e M c 



2500 


-^2500 


= 220 + 55 


O.86+021 x 10 i4 


0.64+°;}f x 10 13 


M-T relation (2) 




(2.4 + 0.6) x 10 14 




Re 


, = 300 


1.38+° 33 x 10 14 


1.09^0 x 10 13 


Y x - M relation (3a) 


495 ± 45 


1.93!°^ x 10 14 


500* 


R500 


= 510 +55 


2.07+°;™ x 10 14 


2.13!°;™ x 10 13 


Y x - M relation (3b) 


560 +5 ° 


2.8^ x 10 14 


200* 


R200 


= 760 +8 ° 


2.8^;° x !0' 4 


2.9 + _q° 7 x 10 13 


NFW (4) 


540 ± 90 


(2.6+ 1.8) x 10 14 (3.7 ± 2.0) x 10 14 



fects on the extrapolated masses. We prefer not to consider this 
value in the following discussion. 

2.5. Systematics in the total mass measurements 

The mass measurements at R^m and R200 have been obtained un- 
der the assumption of hydrostatic equilibrium as an extrapolation 
of the observed profile beyond the largest radius where X-ray 
emission is detected (300 kpc), as well as under some reason- 
able assumptions on the slope of the ICM density and tempera- 
ture profiles beyond this radius. Our assumptions are: i) the slope 
of the density profile remains constant, ~ 2.2 beyond 300 kpc; ii) 
the temperature profile declines as oc r~° 24 outside the extraction 
radius as observed in local clusters in the (0.1 - 0.6) x r^o range 
(Leccardi & Molendi 2008). In order to evaluate possible sys- 
tematics associated to the extrapolation of the total mass to R^m 
and R200, we present the mass estimates adopting three different 
methods. 

The first estimate (method 2) is based on the redshift- 
dependent scaling relations calibrated on local clusters and pre- 
sented in Vikhlinin et al. (2009). From the empirical relation de- 



scribed in their Table 3, we find M 500 = (2.4 + 0.6) x 1O 14 M . 
As a third method, we use the integrated Compton parameter 
Yx = TxxM/cm, which is considered a robust mass proxy within 
7?50o, as shown by numerical simulations (e.g., Kravtsov et al. 
2006). The observed value for CX01415, assuming r^m as ob- 
tained from method 1, is Y x = 1-24+°-^ x 10 14 keV M . 

If we use the Yx-M empirical relation taken from Table 3 
of Vikhlinin et al. (2009), we obtain fl 50 o = (495 + 45) kpc 
and M500 = 1. 93^ 43 x 1O 14 M . Using the calibration based on 
numerical simulations obtained in Fabjan et al. (2011), we find 
M500 = x 1O 14 M , for R 500 = 560^° kpc. Both values 

are obtained iteratively in order to compute consistently all the 
quantities at rsoo- The discrepancy of a factor -1.5 between the 
numerical simulations and the phenomenological Ms^-Yx rela- 
tion can be partially explained with some nonthermal pressure 
support not accounted for in the hydrostatic mass estimates (see 
Kravtsov et al. 2006). 

We also extrapolate the mass according to the Navarro, Frenk 
& White profile (NFW, Navarro et al. 1996, method 4), after re- 
quiring a normalization at 300 kpc to the total mass value ac- 
tually measured from the data. Still, the extrapolation depends 
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Table 3: Summary of the optical/NIR/lR imaging datasets. 
BVRiz are from Subaru- Sup rime, JK S from Subaru-Mo/rci, and 
3.6pm from Spitzer-IRAC. 



Band 


Exposure 


Date obs. 


Image quality 




(s) 




(") 


B 


6x360 


2006-05-30 


1.10 


V 


6x360 


2003-04-05 


0.66 


R 


6x480 


2003-04-04 


0.74 


i 


3x300 


2001-06-27 


0.70 


z 


9x180 


2001-06-27 


0.62 


J 


2160 


2009-04-19 


0.62 


K, 


1440 


2009-04-19 


0.56 


3.6/jm 


97 


2006-02-05 


1.7 



on the unknown concentration parameter. However, such de- 
pendence is not strong and can be included in the extrapolated 
value once we restrain the concentration in the plausible range 
Cnfw = 4.0 ± 0.5, as found in simulations for the wide mass 
range 5 x 10°M o < M < 5 x 10 14 M o (Gao et al. 2008). We find 
M 500 = 2.6^3 x 1O 14 M and M 200 = 3.7+2-0 x 10 i4 Me _ Here 
the much larger error bars include both the statistical error on 
the measured mass at R = 300 kpc and the uncertainty on the 
concentration parameter. 

The systematic uncertainty due to the different mass mea- 
surement methods can be appreciated in Table 2. We find that 
the systematics error on M500 is about 20%, slightly lower than 
the lcr confidence level associated to the statistical uncertainty, 
which is ~ 30%. We note that mass estimates based on a beta- 
profile fitting and the assumption of the hydrostatic equilibrium 
can give values ~ 20% lower around rsoo (and even larger for 
larger radii) due to the violation of hydrostatic equilibrim and the 
presence of significant bulk motions in the ICM. This has been 
shown in numerical simulations (see Bartelmann & Steinmetz 
1996; Rasia et al. 2004; Borgani et al. 2004). Observational cal- 
ibrations also show that hydrostatic masses underestimate weak 
lensing masses by 10% on average ar rsoo (see Hoekstra 2007; 
Mahdavi et al. 2012). We do not find evident bias when compar- 
ing our hydrostatic mass (method 1) with values obtained with 
empirical calibrations (methods 2, 3), also because the statisti- 
cal errors are larger than the possible effects of residual bulk 
motions. Therefore, all the mass measurements we obtained are 
consistent with each other within the errors. Our conclusion is 
that it is feasible to measure M500 for high-z, massive clusters 
using deep Chandra exposures, although clusters with fluxes be- 
low 10~ 14 erg s _1 ctrT 2 , like CX01415, suffer a nonnegligible 
systematic uncertainty due to some unavoidable extrapolation. 



3. Optical/IR analysis 

The field of CX01415 has wide optical photometric cover- 
age with Subaru-SMpr/me (BVRiz), in addition to near-IR JK f 
with Subaru -Moires and Spitzer-IRAC 3.6pm. A summary of 
the archival photometric data used in our analysis is provided 
in Table 3. Standard reduction procedures using IRAF routines 
were used to produce the optical/NIR science images retrieved 
from SMOKA (Baba et al. 2002), with an average seeing of 
0.68" in the optical bands, and 0.60" in the near-IR. Finally, we 
used level 2 Spitzer-IRAC 3.6 pm observations from the archive, 
which are higher level reduced (pbed) products. 



3. 1 . Color-magnitude diagram 

Relaxed galaxy clusters are characterized by a well-defined se- 
quence of red, passively evolving galaxies (see Gladders & Yee 
2005). This relation has been shown to remain fairly constant up 
to z = 1.4 (Mei et al. 2009; Strazzullo etal. 2010; Rosati et al. 

2009) . However, recent studies at z > 1 .4 indicate the appearance 
of star formation in the central brightest galaxies (Hilton et al. 
2010; Fassbender et al. 2011b; Santos etal. 2011), with an in- 
creasing scatter in the red sequence. 

In this section, we investigate the (z' - 3.6pm) color- 
magnitude relation (CMR) of CX01415. This optical/IR diag- 
nostic has proven to be an accurate tool to study the galaxy pop- 
ulation of z > 1 clusters (e.g., Muzzin et al. 2009; Brodwin et al. 

2010) . The BRz color image reveals a population of red galaxies 
centered on the X-ray contours, as shown in Figure 5. We use 
SExtractor (Berlin & Arnouts 1996) to perform source detection 
and galaxy photometry. We measure aperture magnitudes cor- 
rected to 1" in the z-band and apertures of 1.8" in the 3.6pm 
band. The total 3.6pm magnitude is given by the SExtractor 
parameter MAGAUTO. An aperture correction corresponding to 
2.5 x log(1.112) was applied to the IRAC band to match the 
optical magnitudes 3 . 

The z' - 3.6pm CMR, which includes only objects within a 
cluster-centric distance of 1', is presented in Figure 6, left-hand 
panel. The galaxy coincident with the X-ray peak (indicated by 
a magenta star) has a total 3.6pm magnitude of 16.9 and a color 
equal to 5.7. Contrary to what is commonly observed at lower 
redshifts, this galaxy is not the brightest cluster galaxy. A se- 
quence of red galaxies clearly stands out at 5 < z' - 3.6pm 
< 6.5. For comparison we also plot the red sequence of 21 
out of 25 spectroscopically confirmed members of WARPJ1415 
(Huang et al. 2009), which has a median color of 4.4 mag. If 
we further restrict the color cut to 5.5 < z' - 3.6pm < 6.3 and 
the total magnitude to below 18 mag, we obtain a more strin- 
gent locus for the red sequence, described by a linear fit with 
a slope of -0.25. The mean (z' - 3.6pm) color computed for 
the enclosed 12 galaxies is then 5.8 mag, with a spread of 0.21 
mag. This already shows that CX01415 is at considerably higher 
redshift than 1 . The cluster red sequence is well approximated 
by a Bruzual & Chariot (2003) single starburst model with pas- 
sively evolving galaxies with an age of 2 Gyr. These models, 
which accurately reproduce the red-sequence of WARPJ1415 
at z — 1, predict a redshift of ~ 1.5 for CX01415 for a color 
{z' - 3.6pm) = 5.8. This simple redshift estimate is improved in 
the next section by using all eight photometric bands. 

3.2. Photometric redshifts 

Using our eight-band multicolor catalog, we computed photo- 
metric redshifts for the candidate red galaxies with the photo-z 
code zebra (Feldmann et al. 2006) by comparing the observed 
magnitudes with those expected from template spectral energy 
distributions (SED) by Bruzual & Chariot (2003) through a stan- 
dard SED fitting procedure. Our results for CX01415 were cal- 
ibrated against WARPJ1415 at z — 1.03, for which we obtained 
an accurate mean photo-z of (z p hot) - 0.98, with an rms disper- 
sion ~ 0.08 based on the colors of 17 spectroscopic members. 
The comparison of spectroscopic and photometric redshift for 
WARPJ1415 members shows that the photo-z accuracy on the 
average redshift is A<z) = 0.05. Of the 21 CX01415 galaxy can- 



3 http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/iracinstrument 
handbook/28/_Toc296497400 
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Fig. 6: Le/f: (z - 3.6/mi) color-magnitude relation of CX01415. Only objects within 1' distance from the X-ray center are shown 
(black circles). The 21 red cluster candidates (red circles) are selected on the basis of the color cut 5 < z - 3.6/im < 6.5. For compar- 
ison we also plot 19 confirmed members of WARPJ1415 at z — 1 (blue triangles). The magenta star indicates the galaxy coincident 
with the X-ray peak. Right Selection of red, high-z cluster galaxy candidates (white circles), with X-ray contours displayed in green. 
The brightest candidate is marked by a yellow circle. 



didates selected in the CMR, four are outside the FoV of the 
MOIRCS data. A visual inspection of these galaxies forced us to 
discard another four galaxies because of their faint magnitudes 
and blending with nearby objects, which may result in unreliable 
photometry. The histogram of the photometric redshifts is shown 
in Figure 7. Based on the SED fit of 13 red candidates, we find 
for CX01415 (Zphot) - 1-52 ± 0.06, where the 1 <x error is esti- 
mated by extrapolating to z — 1 .5 the photo-z accuracy measured 
at z = 1, assuming a dependence oc (1 + z). This result is in good 
agreement with the zx obtained from the X-ray spectral analysis. 

At present, we do not yet have the spectral confirmation of 
the redshift of CX01415. We recently analyzed long-slit spectra 
of two galaxies (including the central one) obtained with OSIRIS 



(Optical System for Imaging and low-intermediate Resolution 
Integrated Spectroscopy, Cepaetal. 2000) on the 10m facility 
Grantecan at La Palma (Spain). Unfortunately, we collected only 
one out of seven hours originally requested, and the quality of 
the spectra was not sufficient to derive reliable redshifts. More 
efforts to achieve a spectroscopic measure of the redshift of 
CX01415 are ongoing. 

Finally, we note that an in-depth morphological study of the 
cluster galaxies would be extremely interesting. Robust morpho- 
logical classification of cluster galaxy candidates would require 
HST/WFC3 imaging. Our Subaru NIR imaging in very good 
seeing data suggests that about half of the selected red sequence 
galaxies show a distorted late-type morphology, which is not 
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Fig. 7: Histogram of photometric redshifts computed in the clus- 
ter field. We show the photometric redshift distribution of 17 
spectroscopic members of WARPJ1415 at z = 1.03 (blue dashed 
line) and the photometric redshift distribution of 13 candidate 
members selected from the red sequence of CX01415 (red solid 
line). The thinner histogram (solid blue line) marks the spectro- 
scopic redshift distribution of WARPJ1415, while the cyan stripe 
corresponds to the X-ray redshift zx = 1.46 ± 0.025. 

uncommon in clusters at these redshifts (see Hilton et al. 2010; 
Fassbender et al. 201 lb; Santos et al. 201 1). 

4. Discussion 

CX01415 was originally discovered by visual inspection in the 
original 80 ks ACIS-I exposure. 4 In order to estimate the effec- 
tive sky area within which it was found, we carried out a sim- 
ple experiment based on an inspection of the same image at dif- 
ferent exposure times. This experiment allowed us to determine 
the minimum exposure time necessary for a mere detection of 
the extended emission. By trimming the exposure time of the 
ACIS-S and ACIS-I image of CX01415, we found that 20 net 
counts can be considered a firm limit for detection as an extended 
source via visual inspection. This implies that CX01415 could 
have been detected in all the ACIS-I observations with more than 
15 ks and all the ACIS-S observations with more than 12 ks. 
The most conservative choice is to compute the solid angle cov- 
ered by all the Chandra ACIS-I and ACIS-S extragalactic fields 
where CX01415 may have been detected. We browsed the entire 
Chandra archive up to Cycle 12, combined multiple exposures 
of the same region of the sky, and found that a total of about 70 
deg 2 had been observed at the required depth. This must be con- 
sidered as an upper limit, since we do not take into account the 
solid angle lost due to the presence of bright extended targets in 
several extragalactic Chandra fields. 

We can use this extremely conservative choice of sky cov- 
erage to estimate the maximum search volume from where 

4 This also suggests that a "crowd-sourcing" project of the kind of 
Galaxy Zoo (see Lintott et al. 2011) may be very useful if applied to 
Chandra X-ray images. 



CX01415 has been drawn. Using this information in combi- 
nation with the cluster mass function of Tinker et al. (2008) 
with WMAP7 parameters (Komatsu et al. 201 1) 5 , we find that 
CX01415 by itself, after considering the uncertainty in its virial 
mass and accounting for the corresponding Malmquist bias 
(Mortonson et al. 2011), does not create any tension with the 
reference ACDM model. In order to compare our data with 
the theoretical predictions, we computed the mass of CX01415 
for a given overdensity with respect to the background density 
(as opposed to the critical density, which is more commonly 
used) according to the masses defined in Tinker et al. (2008). 
Therefore we measure the mass of CX01415 within an over- 
density of A = 600 with respect to the background density, find- 
ing M 6 oo = 2.0^ x 1O 14 M . Finally, we find that six clusters 
with the mass of CX01415 are expected at z > 1.46 in 70 deg 2 . 
If we consider a 20% higher mass as expected in the presence 
of a bias due to the violation of the hydrostatic equilibrium (as 
discussed in Section 2.5), we find that only three clusters are ex- 
pected at z > 1.46 in 70 deg 2 . Therefore, even if CX01415 is 
not an extreme cluster, it is clear that a deep investigation of the 
entire Chandra archive in search of faint extended sources can 
provide a sizeable number of clusters with masses ~ 1O 14 M . 
Such a sample can provide tight constraints on the cosmolog- 
ical parameters and possibly some discrepancy with respect to 
the reference ACDM. On the other hand, to achieve the same 
accuracy on the mass on the entire sample, a large amount of 
Chandra time would be required. 

The mass determination of CX01415 is also difficult as it ap- 
pears to be less luminous than other distant massive galaxy clus- 
ters. Its total 6 soft-band luminosity is estimated to be Lx[0.5 - 
2.0] = (0.93 + 0.10) x 10 44 erg s _I . Since the measured mass of 
WARPJ1415 at z = 1.03 is comparable to that of CX01415, it is 
interesting to compare their luminosities. We find that the total 
luminosity of WARPJ1415 is L x [0.5 - 2.0] = (4.0 + 0.2) x 10 44 
erg s . If we remove the contribution of the cool core of 
WARPJ1415, using the double beta model fitting of Santos et al. 
(2012), we find a total luminosity of L x [0.5 - 2.0] ~ (3.4 + 
0.2) x 10 44 erg s _1 , which is the correct value to be compared 
with CX01415. Therefore, we find that CX01415 is more than 
three times less luminous with respect to WARP1415. 

It is also interesting to compare these values with those ex- 
pected from the empirical, redshift-dependent M-L relation pre- 
sented in Vikhlinin et al. (2009). We find that WARPJ1415 is 
expected to have Lxexp = (1.6 ± 0.4) x 10 44 erg s _1 , while 
CX01415 Lxexp = (2.0 + 0.5) x 10 44 erg s -1 . The expected 
values do not include an uncertainty of a factor ~ 1.5 due to 
the observed intrinsic scatter in luminosity for a fixed mass. 
Remarkably, WARPJ1415 is scattered above the average L-M 
relation by a factor of 2, while CX01415 is scattered below it 
by about the same factor. If we consider the most massive X- 
ray cluster investigated so far at high z, XMM2235 (Rosati et al. 
2009), we find once again that the measured total luminosity 
L x [0.5-2.0] = (3.0 + 0.15)xl0 44 is lower than the expected one 
by a factor of 2. Finally, the less massive RXJ1252 (Rosati et al. 
2004) at z = 1 .235 is consistent with the average L-M relation. 



5 Specifically, we assumed a flat ACDM model with erg = 0.81, Clo = 
0.275, h = 0.70 for a flat geometry. 

6 The total luminosities are obtained by integrating the surface bright- 
ness distribution up to the nominal virial radius. Although the large ma- 
jority of the emission is concentrated towards the center, the total lumi- 
nosity is estimated to be about a factor 1.4 larger than what was actually 
measured within the extraction region. 
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Fig. 8: Chandra images of four high-redshift massive clusters observed with deep Chandra observations. The four clusters are all 
re-imaged at the same redshift of CX01415 (z = 1.46) and their observations rescaled to a 100 ks ACIS-S observation (including 
CX01415) in order to offer an immediate comparison. Top left: CX01415; top right: WARPJ1415; bottom left: XMM2235; bottom 
right: RXJ1252. Images are smoothed with a 1 arcsec Gaussian kernel and cover 3x3 arcmin. 



The striking difference of the appearance of the four clus- 
ters is shown in Figure 8, where XMM2235, WARPJ1415, and 
RXJ1252 are re-imaged at the same redshift of CX01415 for a 
hypothetical 100 ks exposure with Chandra ACIS-S. Although 
very small, K-correction effects are taken into account, and the 
proper background for a 100 ks ACIS-S observation is simu- 
lated. Therefore, these images can be considered as real data, 
since they were obtained simply by trimming real exposures in 
order to achieve the exact brightness that each cluster would have 
at z = 1-46 for a 100 ks ACIS-S observation 7 . The only artificial 
interventions are the almost negligible resizing (since they are 
all high-z clusters and the angular distance has a weak depen- 
dence on z for z > 1) and the added background in order to 
be consistent with a real ACIS-S observation. Thanks to these 



7 We also note that the profile of CX01415 shown in Figure 8 is 
noisier than that shown in this paper, since it refers only to ~ 1/3 (100 
ks) of the total exposure time. 



properties, the images offer an immediate visual comparison of 
the four clusters as if all of them were at the same redshift of 
z = 1.46. We perform a simple radial profile analysis on the 
four images shown in Figure 9. Clearly, the differences in lu- 
minosity are not only due to the presence of a cool core, but 
also to the overall normalization of the brightness profiles. The 
diversity of this small sample suggests an increasing spread in 
the properties of high-z clusters, possibly related to the transi- 
tion from the formation of protoclusters with little ICM emission 
(z > 2) and the first virialized massive clusters around z ~ 1.5. 
A large fraction of clusters at this redshift may be in the first 
stage of collapse, an occurrence that could explain both the low 
luminosity of CX01415 and its elongated shape. We argue that 
the deep X-ray investigation of very high redshift clusters will 
provide useful insights into the ICM and dynamical properties 
of clusters in this crucial transition epoch. A high-z cluster that 
will be observed with a deep exposure by Chandra in Cycle 14 
is XMMUJ0044 at z = 1.578 (Santos et al. 2011). Another re- 
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Fig. 9: Radial profiles from the four clusters re-imaged at 
Z = 1.5 shown in Figure 8: WARPS 1415 (solid black line); 
XMMUJ2235 (dashed magenta line); RXJ1252 (dotted green 
line); CX01415 (long-dashed red line). 



markable target is the IR-selected cluster IDCS1426 at z — 1.75 
(Stanford et al. 2012), which has a flux five times higher than 
CX01415. Unfortunately, at present this cluster has only a 10 ks 
Chandra ACIS-I exposure and therefore no X-ray mass estimate. 
We note, however, that IDCSJ1426 will be observed for 100 ks 
in Chandra A014 (PI M. Brodwin). 

Finally, as for the chemical properties of the ICM, we no- 
tice that the measured Fe abundance of CX01415 is superso- 
lar, but the 1 <x error is large enough to make it consistent with 
the average values measured at z ~ 1.3 (see Tozzi et al. 2003; 
Balestra et al. 2007). Given the flux of CX01415, an exposure 
time exceeding lMs would be necessary to measure the Fe abun- 
dance with a 20% error and to confirm the Fe excess. This would 
be extremely interesting though, since this value cannnot be 
traced back to the presence of a cool core hosting a brightest 
cluster galaxy (BCG). Both the low value of the csb parameter 
and the absence of a clear BCG in the optical image exclude the 
presence of a strong cool core. If confirmed, the high Fe abun- 
dance value would be another striking difference with respect 
to low-z clusters, where supersolar Fe abundances are always 
associated to cool cores and clear BCG (De Grandi et al. 2004; 
de Grandi & Molendi 2009). In addition, the lookback time of 
9.2 Gyr can put strong constraints on the early phases of ICM 
enrichment by Type la supernovae. This is clearly a case for fu- 
ture X-ray facilities with high effective area and high angular 
resolution 

To summarize, although the number of cluster candidates in 
the range 1 < z < 2 is steadily increasing, mostly thanks to 
IR surveys, very few deep X-ray studies are available. Targets 
drawn from non-X-ray selected samples may turn out to be very 
faint in X-ray, and deep Chandra exposures will remain the 
only means to physically characterize such faint sources. It is 
clear that only a sensitive, survey-dedicated mission can pro- 
vide a major breakthrough in this field. Even though eROSITA 



(Predehl et al. 2010; Cappelluti et al. 2011) will finally provide 
an X-ray all-sky coverage 20 years after ROS AT, it does not have 
the sensitivity or the spectral response needed to find and char- 
acterize clusters at fluxes as low as few x 10~ 14 erg s _1 cm -2 . 
On the other hand, a wide field imager with CCD spectral res- 
olution and an angular resolution as good as 5 arcsec (HEW) 
has been shown to be highly effective in this science case (see 
Murray & WFXT Team 2010; Rosati et al. 2010), and in many 
others concerning the rich, high-z X-ray sky. 



5. Conclusions 

In this paper, we presented the deep X-ray observation of 
the most distant X-ray cluster of galaxies discovered with 
Chandra so far. CX01415 was serendipitously discovered in the 
WARPJ1415 field, and imaged in the deep exposure aimed at 
studying the ICM properties of WARPJ1415 at z = 1.03. Even 
though we do not have optical spectroscopy, our rich multi- 
wavelength observations allowed us to estimate a robust pho- 
tometric redshift, based on the CMR technique and on stan- 
dard SED fitting procedures. Our photometric analysis is fur- 
ther strengthened by the presence of WARPJ1415 in the imag- 
ing fields, which acted as redshift calibrator. The photometric 
redshift of CX01415 turns out to be z p hot = 1.52 ± 0.06. This 
result is confirmed by the detection of the K a Fe emission line 
in the X-ray spectrum, which gives zx - 1.46 ± 0.025. Contrary 
to the few currently known very high-redshift clusters, we are 
able to obtain a robust measure of the ICM temperature and to- 
tal mass within 300 kpc and a reliable extrapolation up to Rsoq- 
Here we summarize our main results: 

- CX01415 is detected with a S /N ~ 11 within a circle with 
a radius of 24 arcsec. The surface brightness can be fitted by 
a beta model; the soft-band flux within r = 24" is equal to 
So.5-2MeV = (6.5 + 0.3) x lO" 15 erg s" 1 cnr 2 ; 

- we identify the Fe K a complex line with a ~ 2.8<x confidence 
level, showing that the measure of redshift through X-ray 
spectral analysis is possible up to the highest redshift where 
X-ray clusters are currently detected; 

- the spectral fit with a mekal model gives kT = 5.8^[ '? keV 
and a Fe abundance Zf e = 1.3+q|Zo; 

- the total mass measured at #2500 = (220 + 55) kpc is M2500 = 
8.6! 2 ' x 10 13 M o . The total mass is measured up to 300 kpc, 
where it amounts to M(r < 300kpc) = 1.38+°|j x 10 14 M o , 
while the ICM mass is M ICM (r < 300kpc) = 1.09+°|° x 
1O 13 M , resulting in a ICM mass fraction of ~ 10%; 

- by fixing the slope of the density profile to the value mea- 
sured at 300 kpc and assuming a mildly decreasing temper- 
ature profile as observed in local clusters, we are able ex- 
trapolate the mass measurement up to R500 = 510^ kpc, 

finding M 50 o = 2.07+° ™ x 1O 14 M ; this value is consistent 
with those obtained using empirical relations between X-ray 
observable and total mass; 

- considering its mass, CX01415 appears to be underlumi- 
nous, particularly if compared with other high-z massive 
clusters; 

- the color-magnitude diagram of cluster member candidates 
shows both a red sequence and a significant fraction of blue, 
irregular-morphology galaxies. A spectroscopic follow up is 
needed to perform an in-depth study of the galaxy popula- 
tion, which is expected to show the signature of the red se- 
quence assembly around z ~ 1.5; 
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- when compared with the expectations for a ACDM universe 
based on the mass function of Tinker et al. (2008), CX01415 
appears to be a typical cluster at z ~ 1.5 for a WMAP 
cosmology (Komatsu et al. 201 1); however, the redshift and 
the total mass of CX01415 place it among the sample of 
massive, distant galaxy clusters which may be used to test 
the standard ACDM model once a complete census of the 
Chandra archive is performed. 

We note that the realm of high-z clusters (z ~ 1 .5 and higher, 
when the lookback time is larger than 9 Gyr) is at the limit of the 
capability of Chandra and well beyond what XMM-Newton and 
eROSITA can reach. At present, only an handful of sources have 
the X-ray coverage sufficient for a detailed analysis. While the 
Chandra and XMM-Newton missions will somewhat increase 
this sample, it is not realistic to foresee a significant increase of 
data on massive high-z clusters with the required high S/N. In 
the near future, this field is likely to be dominated by SZ obser- 
vations and other near-IR large-area surveys, which will unveil 
more and more clusters at z » 1 . However, X-ray observations 
will still be required for a physical characterization of these sys- 
tems and a vast range of astrophysical and cosmological applica- 
tions. Only a wide-field, high angular resolution X-ray mission 
with a large collecting area and good spectral resolution up to 7 
keV seems to be able to match such requirements. 
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